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9 Abstract 

10 Presented is a total- Lagrangian displacement-based non-linear finite-element model of thin-film membranes for sta- 

1 1 tic and dynamic large-displacement analyses. The membrane theory fully accounts for geometric non-linearities. Fully 

12 non-linear static analysis followed by linear modal analysis is performed for an inflated circular cylindrical Kapton 

1 3 membrane tube under different pressures, and for a rectangular membrane under different tension loads at four comers. 

14 Finite-element results show that shell modes dominate the dynamics of the inflated tube when the inflation pressure is 

15 low 1 , and that vibration modes localized along four edges dominate the dynamics of the rectangular membrane. Numer- 

16 ical dynamic characteristics of the two membrane structures were experimentally verified using a Polytec PI PSV-200 

1 7 scanning laser vibrometer and an EAGLE-500 8-camera motion analysis system. 

18 © 2004 Published by Elsevier Ltd. 

19 Keywords: Membrane dynamic characteristics; Non-linear finite elements; Scanning laser vibrometer; Camera-based motion analysis 

system 

21 


22 1 . Introduction 

23 Recently there has been a renewed interest in deployable/inflatable structures for terrestrial investiga- 

24 tions. Moreover, because the cargo space of a launch vehicle is always limited, large space structures must 

25 be designed to be stowed during launch and deployed once on orbit. Hence, instead of using previous elec- 

26 tro-mechanical deployment systems, recent efforts of NASA concentrate on the use of inflatable structures 

27 for space applications (Salama et al., 2000; Jenkins, 2001). For example, NASA conducted the Inflatable 
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28 Antenna Experiment (LAE) in 1996 (Domhiem and Anselmo, 1996). The membrane antenna having an in- 

29 flated diameter of 50 ff with three 93 ft long stmts was transported by the space shuttle Endeavour in a 

30 1-fit x 3 ft x 1 .5 ft container. Also inflatable membrane structures have been used in parabolic antennas, 

31 radiators, solar concentrators, sun shields, habitats, radio-frequency structures, optical communication sys- 

32 terns, radars, lightweight radio-meters, telescopes, etc. Moreover, large balloons are also membrane struc- 

33 tures that have been used for many scientific missions. Advantages of membrane structures include small- 

34 stowed volume, lightweight, low cost, and good thermal and damping properties (Palisoc, 1994). However, 

35 there are difficulties in the design of large scientific membranes (Damle et al., 1997). 

36 Over the last few decades, studying the dynamic behaviors of inflatable membrane structures has proven 

37 to be a challenging job. Many researchers have studied the dynamic characterization of membranes using 

38 numerical methods and, when possible, experimental approaches. Numerical methods such as finite differ- 

39 ence and boundary elements were used by some researchers to compute vibration modes and frequencies of 

40 inflatable dams (Hsieh and Plaut, 1990). The membrane material used in the numerical analysis was as- 

41 sumed inextensible and its weight w’as neglected in the determination of the equilibrium shape. They found 

42 that the membrane’s mass density is of little influence on the computed natural frequencies. Other research- 

43 ers used finite elements and boundary elements to model and compute natural frequencies and mode shapes 

44 of a single-anchor inflatable dam (Mysore and Liapis, 1998). They found that the rigid foundation that an- 

45 chors the dam increases the frequencies whereas the presence of impounded water tends to reduce the fre- 

46 quencies. They noted that the natural frequencies are dependent on the internal pressure as well as the 

47 hydrodynamic pressure of the impounding water. The pressure in an inflatable structure can also play a 

48 critical role in the suppression of vibration (Choura, 1997). This study found that the vibration suppression 

49 of inflatable structures can be accomplished by varying the internal pressure and thus there is no need of 

50 other external actuators for vibration suppression. 

51 Some researchers tested extremely lightweight inflatable structures in a vacuum chamber and in the 

52 ambient atmospheric condition (Slade et al., 2001). They found a lack of correlation between the two cases, 

53 and they explained it to be caused by air damping. Because the coupling of a lightweight membrane and air 

54 is a highly non-linear and localized fluid-structure interaction problem, it is difficult to perform accurate 

55 numerical modeling and simulation of such problems. Hence, testing inflatable structures in vacuum con- 

56 ditions becomes necessary in order to verify numerical predictions. Moreover, because of the size limitation 

57 on actual vacuum chambers, tests in vacuum conditions for large membrane structures are only possible by 

58 using scaled models (Pappa et al., 2001). Johnson and Lienard (2001) obtained the natural frequencies and 

59 mode shapes of a one-tenth scale Next Generation Space Telescope (NGST) using a finite-element model 

60 developed using the cable network method. The difference between predicted and measured natural fre- 

61 quencies ranges from 2% to 27%, and it was noted that predicted mode shapes correlated well for strut- 

62 dominated modes, while membrane-dominated modes showed less correlation. The study of pre-stressed 

63 membranes by Hall et al. (2002) showed that the natural frequencies in air are lower than the ones in vac- 

64 uum because air acts as a non-structural mass. But, the numerical natural frequencies obtained by Kuka- 

65 thasan and Pellegrino (2002) were lower than experimental vacuum ones and the error was attributed to an 

66 inaccurate tension force or Young’s modulus. However, they stated that the error reduced as the tension 

67 force was increased. Experiments also showed that it is difficult to excite global vibration modes of a mem- 

68 brane structure by applying excitations at inflatable components because these components have high local 

69 flexibility, and resonant frequencies may vary with the excitation location (Pappa et al., 2001; Gaspar et al., 

70 2002). Moreover, because the fight weight of membranes, contact sensors cannot be used in testing and 

71 non-contact sensors (e.g., scanning laser vibrometers) need to be used (Gaspar et al., 2002). 

72 In recent years many researchers used commercial finite-element packages to model and analyze non-lin- 

73 ear elastic problems of thin-thickness membrane structures (Wong and Pellegrino, 2003; Kukathasan and 

74 Pellegrino, 2003; Johnston and Parrish, 2003; Jha and Inman, 2003; Holland et al., 2003; Greschik et al., 

75 2003). Because of thin thickness, how to prevent wrinkling becomes the major task in the design of mem- 



L. G. Young el al. / International Journal of Solids and Structures xxx ( 2004 ) xxx-xxx 


3 


76 brane structures (especially those for communication use), and hence membrane wrinkling has been heavily 

77 studied (Wong and Pellegrino, 2003; Kukathasan and Pellegrino, 2003; Jacobson et al., 2004; Su et al., 

78 2003; Sutjahjo et al., 2004). Other important problems studied include air mass effect (Giraudeau et al., 

79 2002; Kukathasan and Pellegrino, 2003), deployment methods and dynamics (Tsunoda et al., 2003), stiff- 

80 ening by internal pressure, follower-force effect by internal pressure (Jha and Inman, 2003), wrinkle-free 

81 design (Sakamoto et al., 2003), crease of membranes, curling of membrane edges caused by residue stresses, 

82 gravity-induced sag and pre-stressing (Jacobson et al., 2004), local buckling due to defects, clefting (Lennon 

83 and Pellegrino, 2000), and experimental techniques for accurate non-contact measurements (Slade et al., 

84 2003; Bales et al., 2003). In wrinkling analysis, the major task is to predict wrinkling regimes, out-of-plane 

85 displacements, wrinkle wavelengths, and dynamic characteristics after wrinkling (Wong and Pellegrino, 

86 2003; Kukathasan and Pellegrino, 2003). Because of thin thickness, such analysis always requires the use 

87 of very small plate or shell elements and hence very long computation time is needed. One can see from 

88 the literature that ABAQUS (2001) is one of the popular finite-element codes used in solving such non-lin- 

89 ear elastic problems (Wong and Pellegrino, 2003; Kukathasan and Pellegrino, 2003; Johnston and Parrish, 

90 2003; Jha and Inman, 2003). However, plate and shell elements in commercial codes usually do not include 

91 effects that are special to membranes, such as air mass effect and pressure-induced follower-force effect. 

92 Moreover, dynamic analysis of post-buckled structures requires a fully non-linear static analysis and then 

93 a linear modal analysis, and it usually requires special attentions in using commercial codes because differ- 

94 ent remedial techniques are used for improving convergence and accuracy, preventing numerical singular- 

95 ity, and/or accounting for large geometric non-linearities. Hence, even the post-buckling analysis results 

96 from popular commercial codes do not always match with experimental data (Kukathasan and Pellegrino, 

97 2003; Jha and Inman, 2003). Furthermore, although some commercial packages can give somewhat reason- 

98 able results, the black-box feeling in using those remedial techniques of commercial packages is always a 

99 hidden pain of researchers. For university professors it is especially painful because it is costly and incon- 

100 venient to use commercial packages for teaching finite-element courses, especially the second finite element 

101 course. Hence, researchers have been improving their in-house finite-element codes by adding new non-lin- 

102 ear elements and solution sequences for analyzing membrane structures (Jacobson et al., 2004; Su et al., 

103 2003; Sutjahjo et al., 2004). 

104 The above review shows that study of membranes is challenging because of the modeling, analysis, and 

105 experimental issues caused by the fight weight, high flexibility, thin thickness, and air mass effect. This work 

106 is to avoid some of these issues by developing a non-linear membrane element by implementing a fully non- 

107 linear membrane theory, to use a scanning laser vibrometer and a motion analysis system for dynamic test- 

108 ing, and to investigate the dynamic characteristics of thin-film membranes subjected to internal pressures 

109 and/or external tension loads. 

110 2. Theoretical background 

111 Fig. 1 shows the three coordinate systems needed for the modeling of initially curved membranes. The 

112 xyz is an orthogonal curvilinear coordinate system with the curvilinear axes x and y being on the unde- 

1 13 formed reference surface of the membrane and the z axis being a rectilinear axis, and the 6/C is an orthog- 

1 14 onal curvilinear coordinate system with the curvilinear axes £ and t) being on the deformed reference surface 

115 and the £ axis being a rectilinear axis. The £ and ij represent the convected configurations of x and y axes. 

1 16 Also, an inertial rectangular coordinate system abc is used for reference purpose in the calculation of initial 

117 curvatures. The i* are unit vectors along the £, t], £ axes, ] k are unit vectors along the x, y, z axes, and i a , i*, 

1 18 and i f are unit vectors along the a, b, c axes. It can be shown that the variations of extension strains e x and 

119 e 2 and shear strain >'6(=>’6i + >’ 62 ) on the reference plane of the membrane are given by (Pai and Young, 

120 2003; Nayfeh and Pai, 2004): 
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Fig. 1 The three coordinate systems for modeling and the undeformed and deformed geometries of a differential reference surface of a 
membrane. 


5ei = T nS/u + T 128/12 + T 138/13 
5<? 2 = ^218/21 + T 228/22 + ^238/23 
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8 /, 


8>’6i = - 


8752 = 


1 36 where 
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cosy 6 (l + e 2 ) 

(1 + e 2 )cosy 62 8y 6 - siny^Se, 4 - siny 62 5 e 2 
( 1 + ei ) cos y 61 4 - ( 1 + e 2 ) cos y 62 

( i 4 - e\ ) cos y 61 5 y 6 4 - sin y 61 6e, - sm y 62 5 e 2 
(1 +e,)cosy 61 + (1 + e 2 ) cos y 62 


5 /, 


<?i = y (1 + u x - vk° 5 + H’^)‘ 4- (y x + uk° 5 +w>(t“ 1 )' 4- (vty — uk° x - vk^Y - 1 
c 2 = y (m.v — vk 4 4- + (l + v y 4 uk 4 4 rvftj ) 4 (w v — uk'Y — /A 2) — i 

>6 = 761 + 762 = sin 1 „r 2 i + T \2T22 4 T 13^23) 


761 


tan' 


(1 4 <? 2 )siny 6 


1 +ci + (1 4 e 2 )cosy 6 ’ 


762 = ^n' 


(1 +e,)siny 6 


4 1 + u x - vk° 5 4 H’k°, 

1 11 = • 


21 


1 -f e\ 

Uy - vk° 4 -f wkl 


62 


4 V, + uk] 4 wk° 6i 

Tu ~ \+7\ ■ 

1 4 lly 4 uk^ 4 U’^2 


1 +e 2 + (1 +e,)cosy 6 


T22 = 


1 4 1 4 €j 

5 /n — 5(1 4 u x — vk° 5 4 wk 0 { ) — bu x — k° 5 bv 4 k^bw 
5 / 12 = 5 (t\ 4 uk° s 4 w'^j) = bv x 4 k° 5 bu 4 k° tl bw 
6/13 = 8 (w x — uk\ — nAr^,) = 5 vt’ x — k\bu — k\ x bv 
5/2i = 8(1/2- — vk° 4 4 w^g 2 ) = bu y — k 4 bv 4 k° i2 bw 
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8/23 = 5 (w v — uk° b 2 — n& 2 ) = 5 w v — k° 62 bu - k 2 bv 
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1 + e\ 

tVy - uk° 62 - vk° 2 
1 + e 2 


(1) 

( 2 ) 


( 3 a) 

( 3 b) 

( 3 ct 


( 4 ) 
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140 Here k° are initial curvatures, u, v, w are displacement components of an arbitrary point on the reference 

141 plane, and u x = cu/dx, u y = cu/cy, etc. The variations of deformed curvatures k t are given by 


(5a) 



'50, 

f 50, ) 

< 50 2 

► - [*,]! 502 

,50 3 

[ 50 3 J 


[ d °' 1 

502 

- [K 2 ] i 502 > 

l S03 J 

[ 503 J 


(5b) 


146 where 


50i = 


cos y 61 


(F 3 i S/21 + T 32 5t22 + T 33 5t 23 


sm y 62 


50, 


60, 


cosy 6 (l +e t ] 

cos y 62 

cosy 6 (l + e,; 


cosy 6 (l +e 2 ) 

sin)’ 6 i ^ 

cosy 6 (l + e 2 ) ' 3 

2 (5) ' 62_5y 6') + ' 2cosy 6 (l 4- e\) 

— sin y 6 r 2 i ^ 5t 2 i + (f \2 — siny 6 r 22 ^5t 22 + — siny 6 7’ 23 ^ 

2 cos y 6 (l + e 2 ) 


(T 31 5f„ 

4- r 32 5r 12 + r 33 5t 13 ) 

(6) 

(r 3I 5/„ 

+ 7’ 32 5r 12 + T 33 5t, 3 ) 

(7) 

)5r 12 + ( 

F 23 siny 6 T i 3 ^5?i 3 



8h 


( 8 ) 



0 

k} 

—k\ 


0 
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—k(,i 

= 

—ks 

0 
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\K,] = 
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—ki 

0 
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k(,i 

0 


_ k(, 2 
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u 


i A U 


T 3 j = (T i2 T 22 — T i 3 7’ 22 )//Joi T 22 = (T \- } T 2 \ — T nT 2i )/Ro, T 33 - (T \\T 22 - T 
Ro = [(f Ufa - f n f 22 ) 2 + (7*, 3 7* 2I T n? 23 ) 2 + (f n f a - f l2 T 2 i) 2 ] W = |cosy 6 i 

155 Using the polar decomposition Jaumann strains B mn can be shown to be (Pai and Palazotto, 1995) 

1 / du . Su . 

Bmn ~ X I X * In 4* "X * lm 

2 \c.r_ ox„ 


(y) 


( 10 ) 


1 59 where u is the local displacement vector relative to the deformed coordinate system c^C, and x, = x, x 2 = y, 

160 and x 3 = z. The fully non-linear strain-displacement relations can be derived by using the concept of local 

161 relative displacements and Fig. 1 to be (Pai and Palazotto, 1995) 

{B} = W ( 1J ) 

1 65 where the change of strains through the thin thickness of a membrane is neglected and 
{£} = {B\\,B22, 2# 12 } t 

{iA} = {(1 + e,)cosy 61 - 1,(1 +e 2 ) cos y 62 - 1, (1 4- e,) sin y 61 + (1 + e 2 ) siny 62 } T 

1 69 For a membrane consisting of orthotropic layers, one can obtain the transformed material stiffness matrix 

170 [g (,) ] for the rth lamina from its principal material stiffness matrix [Q (f) ] and its ply angle (measured with 


( 12 ) 
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171 respect to the axis x) by using tensor transformations and assuming that Jaumann stress = 0. The 

1 72 stress-strain relation for the rth lamina is given by 






w 


176 where 


{>>} 


/ r(i) r( 0 T ( ') 
Win J 22 12 


} T , [e (,) ] 



SB 

o (i) 1 
*-16 

SB 

J& 

0 {i) 
y 26 

.gS« 

We 

«-66 J 


(13) 


(14) 


1 79 3. Finite-element formulation 


180 To study the dynamic characteristics of a tensioned membrane, we first perform fully non-linear static 

181 analysis of the membrane under static loading and then perform a linear modal analysis to obtain natural 

1 82 frequencies and mode shapes with respect to the statically deformed configuration. To derive finite-element 

1 83 equations for fully non-linear static analysis we use the principle of virtual work, which states that ( Wash- 

184 izu, 1982) 

5n = 51T nc (15) 

188 where FI is the elastic energy and lV nc is the non-conservative energy due to external loads. 

189 3.1. Elastic energy 


190 Because the elastic energy' If is due to relative displacements among material particles, we have 


sn = L h - i - )5 Gl d '' '■) + <r> ■ u6 (S dQ i -) + If> ‘■ ,s (M dx> ' *■). 


(16) 


1 94 where f. is the force vector acting on the deformed surface of the undeformed area d.v,„d.\\, (i m n) of 

195 an undeformed infinitesimal cube chqch^ch^ (=dF u ). Using the polar decomposition, Jaumann stresses can 

196 be proved to be (Pai and Palazotto, 1995) 


Jn - 


f, i. 


2Ji 


f ! ‘2 , f 2 *1 , . r 

+ “ ; — — 2J 21 ’21 ~r d 12 


dx 2 cbc3 ’ ~~ lz dr 2 dx3 dxidx3 

200 Using (10), (13) and (17) in (16) one can show that 

r 2 i + 1 


5fl = ^ J J +7^5522 +2J? 2 bB n )<iAdz = J2 J J {^[if’jWd^dz 


(17) 


(18) 


204 where A denotes the undeformed area of the reference surface, N is the total number of layers, and z, and 

205 Zj + 1 indicate the locations of the bottom and upper surfaces of the ith layer. Substituting (11) into (18) 

206 yields 

6FI = J^{b\J/} T [<P]{\l/}dA (19) 

210 where [$] is a symmetric 3x3 matrix given by 
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w=e rp y 

i=i jz < 

213 It follows from (12) and (l)-(4) that 

m = 

217 where 

{ U} = {w, u x , u y , v , v x , v v , w, w x , w v } T 

UJ _ 

u dUj 

223 The non-zero elements of [¥] are: 

V[2 = Cll/Oh ^15 = Cu/Co, V|g = C13/C0, 
V13 = C14/C0, V,8 = C 15 /C 0 , Vi9=C]6/Co, 

V„ = (C,2*S - C, 3 *? + C 15 *J - C\(,k° 62 )/Co 
V 14 = — (CiiA'j + Ci 3 ^, + C 14^4 + C 1^2) /Co 
V17 = (Cn*° + Cn^i 4 - Ci 4^°2 + C, 5 ^)/Co 

V22 = C21 /Co, V25 — Ch/Cqi V28 — C23/C0, 
V23 = C24/C0, ^26 = C25/ Co, V29 = C 2 f,/Co: 

V 2 , = (C22*? - C23*? + C25A:" - C 2 i k%)/C, 

V 24 = — (C 2 ,Jfc? + < 723*61 + + C 26 *S)/Co 

V27 = (C 2 \k\ + < 7*22 *6 1 + ^ 24*62 + C 2 sk\)/C$ 


T 32 = 2 C 31 /C 0 , V 35 = 2C32/C0, T 38 = 2 C 33 /C 0 , 

T33 = 2C34/C0, V 36 = 2C35/C0, V39 = 2 C 3 6 /Co, 
v 3 , - 2 (C 32 ^5 - C 33 *? + Css*; - C 364 )/Co 

V 34 = — 2(C31^5 T ^33*61 T ^ 34^4 + C36*5)/Co 

V 3 7 = 2(C 3 ,*? + C 32 *61 + ^34^2 + C 35 ^)/Co 

230 where 

C 0 = (1 +e,)cosy 5 , + (1 + e 2 )cosy 62 


Cni = 


T 2 1 -siay 6 T u 


T 22 ~ siny 6 7~i2 
COS >’ 6 


f 23 - sin y 6 Ti 3 
cosy 6 



( 20 ) 

( 21 ) 

(22a) 

(22b) 

(23) 

(24) 

(25) 
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( 27 ) 


(28) 


C u = (1 +e ] + (1 + e 2 )cosy 62 cosy 61 )r„ - (1 + e 2 )siny 61 cosy 62 C 0 i 
C \2 = (1 + c, + (1 + e 2 ) cosy 62 cos 76i)^i 2 ~ (1 + c 2 ) sin y 6] cos y 62 Co 2 
Ci 3 = (1 +e, + (1 + e 2 ) cos y 62 cos y 6 i)T - (1 + e 2 ) siny 61 cosy 62 C 03 
Cm = -(1 + e,)siny 61 siny 62 f 2 , - (1 + e,)siny 6 , cosy 62 Co 4 
Ci 5 = -(1 +e,)siny 61 siny^f^ - (1 + e,)siny 61 cosy 62 C 05 
Ci 6 = -(1 + e, ) sin y 6 , siny 62 f 23 - ( 1 + e, ) sin y 6 , cosy 62 Co6 

C 2 , = -(1 + e 2 ) sin y 6 , siny 62 7\, - (1 4- e 2 )siny 62 cosy 6 ,C m 
Cn = -(1 + e 2 ) siny 61 siny 62 r I2 - (1 + e 2 )siny 62 cosy 6 iC 02 
C 23 = -(1 4- e 2 ) siny 6 i sin y 62 7" ]3 - (1 + e 2 ) sin y 62 cos y 6 , C 03 
C 24 = (1 +e 2 + (1 4- e, ) cos y 62 cos y 61 ) T 2 , - (1 4- e x ) sin y 62 cos y 6 , C w 
C 25 = (1 +e 2 + (l + ei) cosy 62 cosy 6 ,)7’ 22 - (1 +ei)siny 62 cosy 6 ,Co 5 
C 26 = (1 +e 2 + (1 + e,)cosy 62 cosy 61 )r 23 - (1 4- e, ) sin y 62 cos y 61 C 06 

C 3 , = (1 + e 2 ) sin y 6 , cosy 62 f „ 4- (1 + e 2 ) cosy 6 , cosy 62 C 0 i 
C 32 = (1 + e 2 ) sin y 61 cos y 62 T i2 + (1 + e 2 ) cosy 6 j cosy 62 C 02 
C 33 = (1 + e 2 ) siny 6 i cos y 62 ?, 3 4- (1 4- e 2 )cosy 61 cosy 62 C 03 
Cm = (1 +ei)cosy 61 siny 62 ? 2 i 4- (1 + e,)cosy 61 cosy 62 Co 4 
C 35 = (1 + ei) cos y 61 siny 62 r 22 + (1 +e,)cosy 6 , cosy 62 C 05 
C 36 = (1 +e,)cosy 6l siny 62 f 23 + (1 + ej)cosy 61 cosy 62 Co6 

239 The way the components of { U] are approximated defines a specific finite element. Using the finite-element 

240 discretization scheme, one can discretize the displacements as 

{u, v, w}' = [Af]{<?^} (30) 

244 where {q J ^\ is the displacement vector of the yth element and [AH is a matrix of two-dimensional shape func- 

245 tions. For example, if it is an eight-node serendipity membrane element, we have 

{? W } = {{? ( * ) } T ,{? (i+1) } T ,-.-,{? (t+7) } T } T (31) 

248 where {q tk> } is the displacement vector of the A'th node and is given by 

{ 9 ( *>} = {«(*),»<*), w (i >} T (32) 

251 Substituting (30) into (22a) yields 

{C/} = [Z)]{^}, [D] = [0][Af] (33) 

255 where [3] is a 9 x 3 matrix consisting of differential operators and [D] is a 9 x 24 matrix. Substituting (21) 

256 and (33) into (19) yields 


(29) 


5n = E [ = {§ q ) 

7=1 J ^ ] 7=1 


'[*]{?] 


( 34 ) 
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260 


264 

265 

266 

267 

268 


271 

272 

273 


276 

277 

278 


282 


286 


290 


294 


298 


301 


where 


l^K^ 1 } = [ P]W[*]W<L4 (35) 

N c is the total number of elements, is the area of the yth element, [l6 ^ ] is the stiffness matrix of the /th 
element, [AT] is the structural stiffness matrix, and { q } is the structural displacement vector. We note that 
[F^] and [K] are not explicitly written. 

Because the structural stiffness matrix is a non-linear function of displacements, the governing equations 
need to be solved by an incremental/iterative method. To derive incremental equations, we let 

{ 9 W} = { 9 0 } + {A 9 M}, {£7} = {t/°} + {Af/} (36) 

where {q 0 } denotes the equilibrium solution and the incremental displacement vector. 

Because {J} = \Q\{ip} and initial strains (or stresses) will be considered in the formulation, we will replace 
{t j/} with 

W-{«M + W (37) 

where {i jj 0 ] represent known initial strains in the initial configuration described by the xyz coordinate sys- 
tem, and {i j/} represent unknown additional strains caused by additional loads. Hence, we obtain the first- 
order expansions of {i//} and [T*] as 

m = Wo) + (O + M{A£/} (38) 

and 


[¥] = pF 0 ] + [H] 

where the entry 3j of [H] is given by 

3 2| A, 


dUjdU k 


A U k 


(39) 


(40) 


Then. w r e use (38) and (39) to expand (35) into a Taylor series and neglect higher-order terms to obtain 
{K^]{q^} = f [[£>3 T ['P°] T [<f>]({^ 0 } + {ib 0 }) + [£>] t ['F 0 j T [^j ['F 0 ] { A L/> + + {^°})W 

(41) 


JAW 


Using (40) and direct expansion, one can prove that 

[S] T [*]({<M + {>A 0 }) = + <A°)} = (* ta + i ° n ) } 

= {wwM*"- + *> u -} 

= [fl]{AI/} (42) 

where [£>] is a symmetric matrix and its entry Q t] is given by 

+ « + « = 

Hence, substituting (42) and (33) into (41) yields 


(43) 
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l* M lk M } = 


K 


{A?«} + [*:«]{*«} 


where 




K j is the so-called elemental tangent stiffness matrix and is given by 
AT W j = jf ^ [£>] t ( ['P°] T [<P] ['P°] + [0])[£>]<L4 


( 44 ) 


(45) 


and 


[ATW]toW}| { ^ }={ ^ } = j^[D]W!*j(ko} + {^°})dL4 


(46) 


We note that 



is symmetric. Eq. (41) shows that initial strains and additional equilibrium strains play 


the same role in the formulation. In other words, any converged equilibrium state can be treated as a new 
undeformed configuration with initial strains {i^ 0 } + ( iA°} , which is s im ilar to the updated Lagrangian ap- 
proach but no coordinate transformation is needed before updating {\j/ 0 } because Juamann strains are de- 
fined with respect to the deformed coordinate system, as shown in (10). 


3.2. External loads 


The variation of non-conservative energy’ due to external loads is given by 
SfV nc = J (r l 8u -!- r 2 bv -f rjSu^cLl = J { 6u 8v 5w }{ ri r 2 r 3 } T d 14 

= = {5<?} T {*} (47) 

y=l 

where r l3 r 2 , and r 3 are distributed external loads per unit area along the axes x, y and z, (.10) has been used, 
{ is the elemental nodal loading vector, [R] is the structural nodal loading vector, and 

{R ul } = j W] T ( n r 2 r, } T d4 (48) 

JAP' 


3.3. Incremental-iterative solution method 

With the use of a loading parameter 2, the governing equation of a static problem can be obtained from 
(3#^ 5) as 

[*]{?} = (49) 

where { i? r }(={ i?}/x) is a reference load vector. Let 

k) = k}, = k°} + {A?},- = k},-, + {k}, (50a) 

= 2 , = ;.° + aa , = + 5/., (sob) 

where i ^ 2, i is the number of iterations in searching for a converged solution when the load increases from 
2 = a 0 to 2 = 2° + A2„ \q°) denotes the equihbrium solution corresponding to a = 2°, denotes the iter- 
ated solution corresponding to 2 = 2„ and {Aq}, denotes the incremental displacement vector correspond- 
ing to the increment A/., = 2, - 2° of the loading parameter. Substituting (50a,b) into (49), expanding the 
results into a Taylor series, and neglecting higher-order terms, we obtain 
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[^],- 1 {S?}, = SA,R},_ 1 + k},_ 1 (51) 

346 where 

{*},_, = tiR},., - ([*]{<?}),_„ W = {*},_, + 

349 Here [K] is a non-trivial matrix only if parametric loadings exist, [K] ; _ , is the local tangent stiffness matrix, 

350 and the sub-indices indicate the state, e.g., [£],_, s [K} {q]={ql i /=; _ To solve for {5?}, and 6/., we use an 

351 alternative version of Riks’ method (Riks, 1979; Lee, 1992) to solve (51) and the following constraint 

352 equation 

= 0 ( 53 ) 

355 which limits the arc-length increment by confining the current search direction to be perpendicular to the 

356 previous accumulated searching direction (Riks, 1979). 


357 3.4. Modal analysis 


358 After the statically deformed geometry of a membrane is obtained using the iteration method shown 

359 above, the tangent stiffness matrix [K] (see (45)) of the final deformed state is the stiffness matrix for the 

360 modal analysis. The mass matrix [M] needed for modal analysis is derived through the variation of kinetic 

361 energy 5Tby using (30) as 

5 T = — f [ pjubu -f tide + wSw)d4 cb = — / { 5<7 “' } T [/v] T /h j/v] {q J '}dA 

Ja Jz j—\ Ja^ 


N e 




(54) 


j= 1 


364 where p is the mass density, [A f JP \ is the elemental mass matrix, [Af] is the structural mass matrix, and 

^ — f — I f (S 

"* — / 4 ^ j— / ,,e v ) 4 

Jz JaW 

367 The elemental tangent mass matrix is the same as [M^\ because it is a constant matrix. 


368 4. Numerical results 

369 4.1. Inflated circular cylindrical tube 

370 We consider a membrane tube having a diameter D = 3 in. and a length L = 47.5 in. with two ends fixed. 

371 The tube was made by overlapping (by 0.5 in.) and gluing two opposite edges of a rectangular Kapton film 

372 using a 0.001 ,/ thick double-sided Kapton tape with silicone adhesive. According to the manufacturer the 

373 Kapton film has Young’s modulus E= 3.7 x 10 5 psi, mass density/? = 2.7552 slugs/ft 3 , thickness 

374 h = 0.002in., and Poisson’s ratio v = 0.34. The seam is assumed to have a width w = 0.5in. and thickness 

375 h = 0.005 in. In numerical analyses we focus on the effects of the longitudinal joining seam and different 

376 pressures on mode shapes and natural frequencies of a pressurized tube. 

377 Fig. 2 shows the tube modeled using rectangular membrane elements with the same thickness, where ele- 

378 ments are shrunk by 20% in order to show the connectivity of elements. Fig. 3 shows a single point Fre- 
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Fig. 3. A single point FRF of the inflated tube without a joining seam when p = 1 .75 psi. where the excitation is at DOF 825 and the 
response at DOF 819 in Fig. 2. 
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quencv Response Function (FRF) with an excitation at DOF 825 (see Fig. 2. a nodal transverse displace- 
ment tv) and the response at DOF 819 when the tube is inflated with p = 1.75 psi. The FRF is obtained 
using the first 100 modes with modal damping ratios assumed to be 0.02 for each mode. We note that clear 

peaks only show at modes #1 (#2) and #3 (#4) and there are no clear peaks in high-frequency ranges be- 
cause of modal coupling. Fig. 4 shows the first 12 mode shapes and corresponding natural frequencies of 
the tube inflated with p = 1.75 psi. Because the thickness of the tube is uniform, mode shapes appear in 
pairs because of the symmetric geometry. Hence, modes #1, #3, #5, #7, #9, and #1 1 are the same as modes 
#2, #4, #6, #8, #10, and #12, respectively. We notice that only modes #1, #3, and #9 are global bending 
modes and most of other modes are local shell modes. Table 1 shows that reducing the internal pressure 
from p = 1.75psi to p — 0.8125psi reduces the natural frequencies, and it also changes the appearance se- 
quence of mode shapes after the eighth mode. For example, the 1 1th mode becomes the new 9th mode, the 
new 11th mode is a global torsional mode that does not appear in pair, and the third bending mode be- 
comes the new 15th (and 14th) mode. Table 1 also shows that, w’hen p decreases, the natural frequencies 
of global bending modes only reduce a little, but the ones of shell modes (e.g.. modes #5 and #7) reduce 
dramatically. 

Fig. 5 show the tube modeled using a thickness of h — 0.005 in. for elements on the seam and 
h = 0.002 in. for other elements. The FRFs shown in Fig. 6 show again that there are no clear peaks except 
the two peaks in the low-frequency range caused by the bending modes. Because shell modes dominate 
high-frequency responses and their frequencies are close to each other, modal coupling is expected and it 
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Mode # 1,2: 97.7 Hz 



z 


Mode #5, 6: 420.1 Hz 



Mode #3, 4: 245.8 Hz 



Mode # 7. 8: 429.6 Hz 



Mode #9, 10:439.5 Hz Mode# 11. 12: 456.8 Hz 



Fig. 4. The first 12 mode shapes and natural frequencies of the inflated tube without a joining seam and having an inflation pressure 
p = 1 .75 psi. 


398 will be difficult in obtaining mode shapes in experiments. Fig. 7 shows four low-frequency mode shapes and 

399 natural frequencies when p = 1.75 psi. Because of the joining seam, the geometric symmetry is broken and 

400 natural frequencies and co 2 are different from co 2 and oj 4 , respectively. However, the first four modes are 

401 still global bending modes, and shell modes still dominate high-frequency responses. We note that co 2 > u>\ 

402 because mode #2 involves the bending of the joining seam, the ct> 2 is higher than the u> 2 in Fig. 4, and the a>, 

403 is lower than the co, in Fig. 4. Modes #2 and #4 do not cause peaks in Fig. 6 because they are bending on 

404 the xy plane (see Fig. 5). 
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Table 1 

Comparison of natural frequencies when the inflation pressure is reduced from p — 1.75 psi to p = 0.8 1 25 psi and without the joining 
seam 



(Hz) 

o) 3 (Hz) 

o) 5 (Hz) 

o) 7 (Hz) 

0)9 (Hz) 

p= 1.75 psi 

97.7 

245.8 

420.1 

429.6 

439.6 

p = 0.8 1 25 psi 

96.7 

244.4 

287.5 

300.6 

438.1 (o) ls ) 



Fig. 5. The inflated circular cylindrical tube with a longitudinal joining seam. 
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Fig. 6. The single point FRF of DOF 1425 with an excitation at DOF 1419 of the inflated tube with a joining seam when p = 1.75 psi. 


405 4.2. Tensioned rectangular membrane 

406 The 22 in. x 23.25 in. x 0.002 in. Kapton membrane shown in Fig. 8 has the same material properties as 

407 the inflated tube presented in Section 4.1, and the tension forces are aligned along the two diagonal lines. 

408 The tension force at each comer is applied through a 1 in. x 1 in. thin aluminum plate glued to the Kapton 

409 film. In the non-linear static analysis, the center point is fixed. In the linear modal analysis, the four comers 

410 are fixed by the four aluminum plates. Fig. 9 shows some low-frequency mode shapes and natural frequen- 

411 cies when T = 1.5 lbs. We note that most modes are local vibration modes around the edges except a few 
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Mode# 1: 13 .4 Hz Mode # 2: 17.7 Hz Mode # 3: 17.7 Hz 



Mode # 6: 20.0 Hz Mode# 8: 21 .19 Hz Mode # 9: 22.6 Hz 



Mode #10: 24.4 Hz Mode # 1 1 : 24.7 Hz Mode #12: 27.0 Hz 



Fig. 9. Some low-frequency mode shapes and natural frequencies of the tensioned rectangular membrane with T — 1 .5 lbs. 

412 global modes. These local modes are due to the non-uniform tension field over the entire membrane with 

413 the four edges being under small tension forces. For example, modes #2 and #3 are apparently local modes 

414 and modes #1, #6, and #9 are more like global modes. We note that modes #2, #3, #6, and #10 are dom- 

415 inated by vibrations of the two edges parallel to the x axis, and modes #4, #5, #7, and #1 1 are the corre- 

416 sponding local modes dominated by vibrations of the two edges parallel to the y axis. Because the 

417 membrane shape is non-square, the tension force along the two edges parallel to the y axis is higher than 

418 that along the two edges parallel to the x axis and hence the natural frequencies of modes #4, #5, #7, and 

419 #11 are slightly higher than those of modes #2, #3, #6, and #10, respectively Table 2 compares the natural 

420 frequencies when T — 21bs with those when T — 1.5 lbs. We note that increasing the tension lorce by 0.51b 

421 does not change the appearance sequence of the first few mode shapes, but the natural frequencies increase. 

422 Fig. 10a and b show FRFs of a response point at the center (DOF 1263, see Fig. 8) and a point close to the 

423 center of an edge (DOF 525), respectively. The FRFs are obtained using the first 70 modes with modal 

424 damping ratios assumed to be 0.02 for each mode. Fig. 9 show’s that the dynamics around the edge is dom- 

425 inated by local modes and the dynamics around the membrane center is dominated by global modes, and 

426 Fig. 10a, b exactly show this phenomenon. Similar to those observed in the FRFs of the Kapton tube, the 

427 local modes have natural frequencies close to each other and hence modal coupling exists in all high-fre- 

428 quency ranges. 


Table 2 

Comparison of natural frequencies of the rectangular membrane when the tension force is changed from T = 1.51b to T = 21b 
to j (Hz) to 2 (Hz) (i) 3 (Hz) <o 4 (Hz) w 5 (Hz) 


T= 1.51b 13.4 17.7 17.7 18.5 18.6 

T=21b 15.5 20.5 20.5 21.4 21.5 
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Fig. 10. Single point FRFs of the tensioned membrane when T = 1 .5 lbs: (a) measured at DOF 1263 in Fig. 8. and (b) measured at 
DOF 525. 

429 5. Experimental results and comparison 

430 5.1. Inflated circular cylindrical tube 

431 Fig. !! shows the circular cylindrical Kapton tube set-up with 100 circular retro -reflective markers on it 

432 and a Ling Dynamic LDS V408 shaker attached to the plastic part that sealed and supported the left end of 

433 the tube. Because the membrane is transparent, the retro-reflective markers were used to make the meas- 

434 urement using a Polytec PI PSV-200 scanning laser vibrometer possible and to enhance the laser signal. 

435 However, the markers are small («0.06in. 2 ) and mass loading from the markers is negligible. The retro- 

436 reflective markers look big in Fig. 11 because they reflected the flashlight from the camera. The seam is 

437 on the opposite side of the 100 measurement markers. Experiments were performed for inflation pressures 

438 p = 0.8125psi and p = 1 .75 psi. For the FFT acquisition in using the scanning laser vibrometer, a 0-1 kHz 

439 periodic chirp excitation was used with 6400 FFT lines. Fig. 12a shows the averaged FRF of the tube with 



Fig. 11. The inflated circular cylindrical Kapton tube with retro-reflective circular markers. 
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Fig. 12. The FRFs of the inflated tube with an inflation pressure p = 1 75psi: (a) averaged FRF, (b) FRF measured at center, and (c) 
FRF measured near the excitation location. 


440 an inflation pressure p = 1 .75 psi. Fig. 12b and c shows FRFs measured at the center (see Fig. 11, #2) and 

441 near the excitation (#1), respectively. We notice that the FRFs have many small peaks but the dominants 

442 peaks are in the low-frequency range, as expected. Moreover, because the membrane is so flexible, it was 

443 difficult to excite high-frequency modes and hence the peaks in high-frequency ranges have small ampli- 

444 tudes. The first six Operational Deflection Shapes (ODSs) corresponding to the peaks in Fig. 12 are shown 

445 in Fig. 13. The ODSs #1 and #2 in Fig. 13 correspond to and agree well with modes #2 and #4 in Fig. 7. 

446 However, the experimental frequencies are lower than the numerical ones by 11.5 Hz (11.6%) and 20.6 Hz 

447 (8.3%), respectively. The discrepancy could be the effect of air mass or due to incorrect material properties 

448 used in the numerical computation. Since Young’s modulus of Kapton films is known to be a function of 

449 strains, temperature, and time, the value of E used in the finite-element analysis may not be accurate. Air 

450 mass may significantly change the dynamic characteristics of membrane structures and needs to be ac- 

451 counted for in the numerical analysis in order to have results comparable to experimental ones (Kukatha- 

452 san and Pellegrino, 2002). 

453 We note that ODS #4 corresponds to a traveling wave and thus is a coupled mode. The ODS corre- 

454 sponding to the peak at 60 Hz was found to be a non-uniform traveling shape, and hence it is not an 

455 ODS at all. It w'as believed to be due to the electric power used to drive the electro-mechanical shaker. 

456 The ODSs in Fig. 13 look like pure bending modes. However, because the measurement was along a line, 

457 it is difficult to know if a measured ODS is a bending mode or a shell mode. 
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Fig. 13. The first six ODSs and natural frequencies of the inflated Kapton tube with a pressure p = 1 .75 psi. 


458 When the internal pressure is reduced from p = 1 .75 psi to p = 0.8125 psi, the frequency of ODS #2 de- 

459 creased from 226.6 Hz (see Fig. 13) to 211.4Hz. 

460 5.2. Tensioned rectangular membrane 

461 Fig. 14 shows the test set-up of the 22m. x 23.25 in. x 0.002in. rectangular membrane. Experiments were 

462 performed with an excitation at a single comer and with a simultaneous excitation at the four corners, 

463 respectively. The number of measurement points is 13 x 16. as shown in Fig. 14. The small circular 

464 retro-reflective markers have a radius of 0.14in., but they look big m the picture because they reflected 

465 the camera flashlight. Because of the local flexibility of membranes, the excitation pomts were positioned 

466 at where the tension cables were connected to the membrane at the four comers. Fig. 15 shows the averaged 

467 FRF and a single-point FRF of the rectangular membrane with T= 1 .5 lbf and an excitation at the lower 

468 left comer. The peak at 60 Hz in Fig. 15a was caused by the 60-Hz electrical power used to drive the shaker. 

469 The many small peaks in Fig. 15b around high-frequency ranges are due to local modes around the edges, 

470 as explained and shown in Figs. 9 and 10. Fig. 16 shows the first four ODSs and frequencies under the single 

471 point excitation. Fig. 17 shows the first four ODSs and frequencies obtained with T — 1.5 lbf and a simul- 

472 taneous excitation at the four comers. In order to simultaneously excite the four comers, the four comers 

473 were supported by a frame and a Z-shape rod was used to connect the frame to the shaker (see Fig. 14). The 

474 obtained FRFs are rougher than those in Fig. 15, which is probably due to the vibration of the rod. 

475 Figs. 16 and 17 show that the excitation location may affect the natural frequencies and mode shapes. 

476 ODSs #1 and #4 in Fig. 16 correspond to ODSs #1 and #3 in Fig. 17, and ODS #3 in Fig. 16 is similar 

477 to ODS #2 in Fig. 17. The ODSs #1, #2, and #3 in Fig. 17 correspond to and agree with the modes #1, #6, 

478 and #9 in Fig. 9. However, the experimental frequencies are much lower than the numerical ones. Since the 

479 rectangular membrane has a large area in contact with the ambient air, air mass significantly lowers the 

480 experimental natural frequencies. Moreover, Fig. 15 show's that there are no clear peaks beyond 20 Hz. 
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Fig. 14. The rectangular Kapton membrane excited at the four comers by a frame with a tension force T = 1.5 lbs at each comer. 



Fig. 15. FRFs of the rectangular membrane with T = 1 .5 lbs when the excitation is at the lower left comer: (a) averaged FRF. and (b) 
single point FRF measured at the center of the bottom edge. 


481 In other words, high-frequency modes are highly coupled, and most experimental ODSs were observed to 

482 be traveling modes due to modal couplings. Hence, it is difficult to obtain experimental high-frequency 
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483 ODSs to be compared with the corresponding numerical mode shapes. The results show that only low-or- 

484 der natural frequencies and ODSs were able to be obtained experimentally because of high structural flex- 

485 ibility, high modal density, severe modal coupling, and heavy air mass. Moreover, experimental natural 

486 frequencies are all lower than numerical ones. This discrepancy is attributed to the influence of heavy 

487 air mass. 

488 A simple way to account for the air mass effect in the vibration of the rectangular membrane is to in- 

489 crease the membrane’s mass density. If the density of the Kapton membrane is increase by a factor of 

490 2.6 to account for air mass when T = 1 . 5 lbf, the numerical natural frequencies of modes #1, #6, and #9 

491 decrease from 13.4, 20.0, and 22.6Hz (see Fig. 9) to 8.3, 12.4, and 14.0Hz, respectively. We note that, 

492 although the 6th and 9th numerical natural frequencies become close to the experimental ones in Fig. 17 

493 (ODSs #2 and #3), the first natural frequency does not match because the air mass effect is not exactly 

494 the same as adding mass to the structure. Hence, advanced analysis by considering the air-membrane inter- 

495 action needs to be performed in order to accurately predict natural frequencies of membranes in air (Kuka- 

496 thasan and Pellegrino, 2003). 

497 Because the membrane thickness is usually so thin, any significant amplitude of harmonic excitation may 

498 cause the vibration amplitude to be larger than the thickness, and non-linear modal coupling, modulation 

499 of several linear modes at an unknown frequency, and/or even chaotic vibration may exist (Nayfeh and Pai, 

500 2004). In vibration testing using a scanning laser vibrometer, if the structural vibration is steady and peri- 

501 odic with a known period T and the recording at each location is controlled by triggering to begin at nT (n 

502 is an integer) after the beginning recording time of the previous measurement point, the actual velocity pro- 

503 file at time t = t k will be the distribution of the measured velocities of all points at nT + t k , where n is dif- 

504 ferent for each point. The ODS corresponding to the velocity profile can be calculated as the velocity profile 

505 divided by Q(—2n/T) only if the vibration is harmonic. Hence a scanning laser vibrometer cannot measure a 

506 transient ODS or even a steady-state ODS having an unknown period. Fortunately a camera-based motion 

507 analysis system can solve these problems because it simultaneously traces all markers and provides truly 

508 Lagrangian descriptions of particle motions. 

509 To check whether the ODSs shown in Figs. 16 and 17 arc real ODSs we also used our new EAGLE-500 

510 real-time motion analysis system shown in Fig. 18 to measure the ODSs of the membrane with T = 1.5 lbf. 

511 In order to have non-localized and non-contact excitations we used a lightweight ISin. x 18.25 in. x 0.09 in. 

512 composite plate set-up to be parallel to the membrane and at 2.5 in. away from the membrane’s center to 

5 1 3 push the surrounding air to excite the membrane. The composite plate was fixed on a Ling Dynamic LDS 

514 V408 shaker, and the excitation amplitude was controlled to be 3.5mm, 4.5mm, and 5.0mm for ODSs #1, 

515 #2, and #3, respectively. Experiments showed that this is a very efficient way of exciting the membrane be- 

516 cause the air mass effect is so significant. The motion analysis system uses 8 high-resolution CMOS (com- 

5 1 7 plementary metal-oxide-semiconductor) cameras to capture pictures of a structure when 8 visible red LED 

518 strobes light up retro-reflective markers on the structure. The cameras and strobes are synchronized to 

519 work at a speed between 0.1 and 2000 FPS (frames per second). For a frame rate between 0.1 and 480 

520 FPS, a full resolution of 1280 x 1024 pixels is used. For a frame rate between 480 and 2000 FPS, a reduced 

521 resolution is used. Using triangulation techniques and the known focal lengths (after calibrations using an 

522 L-frame with 4 markers and a T-wand with 3 markers) of the cameras and the known coordinates of the 

523 bright points (caused by the retro-reflective markers) on the 2D pictures inside the cameras, the EAGLE 

524 real-time software EVaRT 4.2 automatically computes and records the instant 3D coordinates of the center 

525 of each retro-reflective marker that is seen by at least two cameras. Hence, 3D time traces of all makers are 

526 available for performing dynamic animation using stick figures and showing pop-up graphs of displace- 

527 ments, velocities, and accelerations, and they can be output to other programs for further signal processing. 

528 The recording time length is effectively infinite and up to 600 markers can be simultaneously traced due to 

529 the use of large computer memory and a 100 Mbit data upload rate. Because the 3D coordinates of each 

530 marker are checked and calibrated when more than two cameras see the marker, the measurement accuracy 




Fig. 18. A typical set-up of an EAGLE-500 motion analysis system for measuring instant coordinates of many markers on a structure. 
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Fig. 19. The first three ODSs and frequencies of the rectangular membrane with T = 1 .5 lbs and excited by the surrounding air pushed 
by a lightweight composite plate. 


531 is high. For example, the measurement error is far less than 1.0 mm when the measurement volume is 

532 2 x 2 x 2m ? . We note that, although the measurement accuracy of the motion analysis system is lower than 

533 that of the scanning laser vibrometer, the ODSs measured by the motion analysis system are real ODSs 

534 because all points are measured at the same time. Fig. 19 shows the first three ODSs measured by the mo- 

535 tion analysis system, which agree well with those in Fig. 17. However, in order to efficiently excite the anti- 

536 symmetric ODS #2 shown in Fig. 19, the excitation plate was moved to excite only the left half of the mem- 

537 brane, but the turbulent air flow around the left edge disturbed the ODS shape. Moreover, the ODS #4 in 

538 Fig. 17 could not be obtained, which is probably due to the way of excitation or it is a non-linear mode with 

539 internal resonance. This problem requires further studies. 


540 6. Concluding remarks 

541 In this paper, we used Jaumann strains and stresses to derive a total-Lagrangian finite-element model of 

542 membranes. Results from finite-element analyses of an inflated circular cylindrical Kapton tube and a ten- 

543 sioned rectangular Kapton membrane were verified by experiments using a scanning laser vibrometer and a 

544 motion analysis system. Finite-element analyses of the inflated tube showed that shell modes dominate the 
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545 dynamics of the tube, and they are difficult to measure using the scanning laser vibrometer. Finite-element 

546 results and experimental data showed that the dynamics of the rectangular membrane is complex. Although 

547 some low-frequency mode shapes correlate well with the experimental ones, numerical natural frequencies 

548 are far higher than the corresponding experimental ones because of heavy air mass that was not included in 

549 the numerical analysis. Moreover, it is almost impossible to verify high-frequency mode shapes because of 

550 high modal density, severe modal coupling, and heavy air mass. 
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